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Abstract 

Population level response to hypoxia has become an issue of global significance 
because of increased frequency and intensity of hypoxic events worldwide, and 
the potential for global warming to exacerbate hypoxic stress. In this study, we 
sequenced two nuclear intronic regions and a single mitochondrial region 
across seven populations of the African cyprinid, Barbus neumayeri from two 
river drainages in Uganda: the Rwembaita Swamp-Njuguta River System and 
the Dura River. We then examined two indices of population structure, G S t 
and Jost's D, to detect links between oxygen availability and genetic variation 
and to determine if population genetic structure was associated with (i) 
dissolved oxygen regime (hypoxia or normoxia), (ii) geographical distance, or 
(iii) a combination of dissolved oxygen regime and geographical distance. Our 
results indicate that over a large scale (between drainages), geographical 
distance significantly affects the genetic structure of populations. However, 
within a single drainage, dissolved oxygen regime plays a key role in determin- 
ing the genetic structure of populations. Within the Rwembaita-Njuguta system, 
gene flow was high between locations of similar oxygen regimes, but low 
between areas characterized by divergent oxygen regimes. Interestingly, G ST 
analyses appear to yield less realistic measures of population structure than 
Jost's D, suggesting that caution must be taken when interpreting and compar- 
ing the results from different studies. These results support the idea that aquatic 
dissolved oxygen can act as a selective force limiting gene flow among popula- 
tions of aquatic species and therefore should be considered when implementing 
conservation plans and assessing environmental impact of human activities. 
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Introduction 

Both geographical isolation and environmental differences 
can decrease population connectivity and contribute to 
genetic divergence between populations, particularly in 
fragmented landscapes (Hedrick 2000; Manel et al. 2003; 
Richardson 2012). Environmental factors can create barri- 
ers to migration and gene flow, potentially fixing life 
history changes, and exacerbating population fragmenta- 
tion and isolation caused by distance. The role of environ- 



mental factors in genetic isolation, and the possibility that 
these factors may trump, exacerbate, or confound diver- 
gence resulting from geographical separation, has received 
renewed interest with environmental factors changing in 
response to the rise of global temperatures, and the sec- 
ondary effects associated with global warming such as oxy- 
gen limitation of thermal tolerance (Davis and Shaw 2001; 
Gienapp et al. 2008; Portner 2010). For example, Qian and 
Davies (1996) showed that interpopulation genetic differ- 
ences between populations of freshwater leeches (Nephelop- 
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sis obscurd) residing in lakes with different dissolved oxy- 
gen (DO) regimes, were more consistent with oxygen act- 
ing as a selective agent than geographical distance. 
Similarly, at least two populations of the common skate 
(Dipturus batis) have been found to be reproductively iso- 
lated around the British Isles (Griffiths et al. 2010). 
Although both skate populations had the same depth and 
substrate preferences, reproduction and dispersal between 
the groups were limited due to areas of low water tempera- 
ture, suggesting that differences in environmental tempera- 
ture can cause genetic isolation of populations in an 
otherwise continuous environment (Griffiths et al. 2010). 
These are not unique examples, and several other systems 
also suggest that environmental differences can cause 
genetic isolation among populations (Reusch et al. 2001; 
Lecomte and Dodson 2004; Bekkevold et al. 2005; Funk 
et al. 2005; McCairns and Bernatchez 2008). 

Several recent studies have demonstrated strong patterns 
of interdemic variation in fish populations between low- 
and high-oxygen environments, particularly in traits 
related to oxygen uptake such as gill size (Chapman et al. 
1999; Timmerman and Chapman 2004; Chapman 2007; 
Binning et al. 2010; Tobler et al. 2011), and structural ele- 
ments surrounding the gills (Chapman et al. 2008), but 
also brain size (Chapman and Hulen 2001; Chapman et al. 
2008), and body shape (Langerhans et al. 2007; Crispo and 
Chapman 2011). These studies suggest that alternative dis- 
solved oxygen regimes may represent a strong source of 
divergent natural selection that may lead to reduced gene 
flow and, potentially, ecological speciation. However, the 
degree to which genetic differences underlie the phenotypic 
divergence across dissolved oxygen regimes appears to vary 
across species. In their studies of the African cichlid Pseud- 
ocrenilabrus multicolor victoriae, Crispo and Chapman 
(2008) examined several populations from the Nabugabo 
and Mpanga River drainages of Uganda to investigate the 
potential for reduced gene flow between populations diver- 
gent in behavioral, morphological, and physiological traits. 
While the authors were able to detect genetic differences as 
a function of geographical distance using one mitochon- 
drial and 12 microsatellite markers, no correlation between 
genetic divergence and the oxygen regime (hypoxic vs. 
normoxic) was observed, suggesting that phenotypic plas- 
ticity may be facilitating migration between ecologically 
different habitats in this system. Common garden rearing 
experiments in this same cichlid species showed evidence 
for very high levels of plasticity in gill size, and both plastic 
and genetics effects in brain mass (Crispo and Chapman 
2010). A second widespread African fish, the cyprinid Bar- 
bus neumayeri that inhabits the higher order tributaries 
and upstream sectors of the Mpanga River drainage dis- 
plays similar differences in behavioral, morphological, and 
biochemical traits between high- and low-dissolved oxygen 



sites (Chapman and Liem 1995; Olowo and Chapman 
1996; Schaack and Chapman 2003; Martinez et al. 2004). 
However, in contrast to the cichlid example, a combination 
of long-term acclimation and allozyme studies suggest the 
possibility of genetic differentiation between oxygen 
regimes over small spatial scales in this species (Chapman 
et al. 1999; Martinez et al. 2004). In this study, we used 
nuclear introns (noted to be variable between fish popula- 
tions, Quattro and Jones 1999; Quattro et al. 2001) to 
examine whether populations of B. neumayeri are geneti- 
cally structured by (i) dissolved oxygen regime (hypoxia or 
normoxia), (ii) physical barriers such geographical dis- 
tance, and/or (iii) both oxygen regime and geography. 
Nuclear introns accumulate change at a relatively fast rate 
(e.g., compared to coding regions) and can often be easily 
amplified from a species of interest using primers designed 
from the relatively conserved flanking coding regions. 

Valley papyrus swamps in East and Central Africa 
provide an excellent system for exploring the role of 
dissolved oxygen as a divergent selective agent, because the 
dense interior of the papyrus swamps is hypoxic through 
the year (Chapman et al. 1999; Chapman and Hulen 
2001), and these swamps represent patches of hypoxic 
habit in the matrix of the larger, flowing, and generally 
well-oxygenated rivers. In papyrus swamps, conditions of 
low light incidence caused by the height of papyrus plants 
(often up to 5 m) and low water mixture, in combination 
with the high amount of organic decomposition of papyrus 
and other plants, result in DO levels well below normoxic 
conditions (Chapman and Liem 1995; Chapman et al. 
1999). Valley papyrus swamps can significantly reduce river 
valleys DO levels for several kilometres, potentially creating 
a barrier to hypoxia tolerant species. We use the largest 
papyrus swamp-river system in Kibale National Park, 
Uganda (Rwembaita Swamp), as the focus of our study 
because B. neumayeri is widely distributed in the park 
occurring in the dense interior of the papyrus as well as 
small tributary streams and the open waters of the larger 
rivers (Langerhans et al. 2007). Kibale National Park, a 
mid-altitude moist forest in equatorial Uganda, is drained 
by two major river systems (Dura and Mpanga) and associ- 
ated swamp and stream tributaries that produce the mosaic 
of dissolved oxygen conditions in which B. neumayeri is 
very abundant. The substantial variation in DO between 
and among ecosystems in Kibale National Park represents 
a system of multiple, semi-isolated, environmental regions 
in which divergent natural selection may be selecting for 
individuals tolerant to a specific dissolved oxygen regime. 

Barbus neumayeri exhibits strong patterns of pheno- 
typic variation between low- and high-dissolved oxygen 
sites within the rivers of Kibale. Fish from hypoxic swamp 
sites are characterized by larger gills (Chapman et al. 
1999; Langerhans et al. 2007), lower thresholds for 
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aquatic surface respiration (Olowo and Chapman 1996), 
and lower critical oxygen tension (Chapman 2007) than 
conspecifics from nearby well-oxygenated stream and 
river sites. Moreover, individuals from hypoxic areas pres- 
ent a higher glycolytic capacity than fish from normoxic 
areas, even after acclimation and acclimatization periods 
(Martinez et al. 2004, 2011). The nonplastic differences in 
enzyme activity in fish from environments with different 
dissolved oxygen levels suggest that these differences are 
genetically based and that the populations of B. neumay- 
eri residing in divergent oxygen regimes are genetically 
isolated. To test this potential isolation, we quantified 
genetic structure using two different methods, G S t and 
Jost's D. G S t the more established metric, quantifies the 
standardized variance in allele frequencies among popula- 
tions. Jost's D coefficient, a more recent derived metric, is 
a measure of genetic differentiation (Jost 2008; Raeymae- 
kers et al. 2012). Some recent work suggests that G S t ma y 



under-estimate genetic differences and that in some cases, 
Jost's D will give a more accurate picture of the true 
genetic structure of populations (Jost 2008; Whitlock 
2011). Here, we used both measures to test the genetic 
differentiation between B. neumayeri subpopulations 
sampled in the Mpanga and Dura River regions. 

Methods 

Study location description 

In this study, we sampled fish from seven locations in the 
Kibale National Park in Uganda, Africa. Sites were selected 
to maximize variation in dissolved oxygen (DO) and geo- 
graphical distance. The two main rivers Dura River and 
the Mpanga River, both eventually feed into the Lake 
George system (Fig. 1). Basic physico-chemical characteris- 
tics of these sites including water current, water tempera- 



Njuguta River 




ikana Stream 



Figure 1. Map of the Dura and Mpanga drainages of the Lake George region, Kibale National Park, Uganda. Inset in the map of Uganda 
indicates the location of the sampling regions. 
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ture, dissolved oxygen concentrations, maximum depth at 
capture locations, conductivity, pH, and water transpar- 
ency are summarized in Langerhans et al. (2007). Six loca- 
tions were within the Rwembaita Swamp-Njuguta River 
system (which feeds into the Mpanga River), and one loca- 
tion was located on the Dura River, deep within the forest. 
The Dura River location is located more than 100 km 
away on a separate tributary of Lake George, and is geo- 
graphically located south of all of the other locations in the 
Rwembaita Swamp-Njuguta River system (Fig. 1). To dis- 
perse from the Dura River to the Mpanga River, fish would 
have to move downstream through a major papyrus 
swamp, and upstream through several kilometres of fast- 
flowing river, and transverse the Mpanga Falls (approx. 
50 m in height). Therefore, contemporary gene flow 
between the Dura River location and the Rwembaita/Njug- 
uta River locations is very unlikely, and a high amount of 
divergence between those locations was expected. The 
seven locations were characterized by divergent DO levels, 
varying from 1.35 mg 0 2 L _1 to 7.32 mg 0 2 L _1 (Langer- 
hans et al. 2007). Selected locations included two swamp 
water locations, two river locations, and three stream loca- 
tions (Fig. 1). The most western location in the Rwembait- 
a-Njuguta system is the Mikana Stream, which flows 
eastward into the Rwembaita Swamp. Two additional 
small tributary streams that feed into the Rwembaita 
Swamp were sampled: Inlet Stream West with Inlet Stream 
East. A tributary of swamp water leaves the Rwembaita 
Swamp (referred to as Rwembaita/Njuguta Outflow 
Stream), and passes through secondary forest for 80 m 
before meeting the waters of the Njuguta River, the latter 
drains into the Mpanga River. DO levels in the Rwembaita 
Swamp have been recorded since the early 1990s, and DO 
has remained remarkably low over the years. Seasonal vari- 
ation is apparent in the swamp water with higher DO dur- 
ing the rainy periods; however, the swamp water remains 
well below saturation throughout the year (Chapman and 
Liem 1995; Chapman et al., 1999; Langerhans et al. 2007; 
Joyner-Matos and Chapman, 2013). 

Fish sampling and DNA extraction 

Fish were sampled during the rainy season in January 2010 
using baited minnow traps, and body mass and total length 
were measured in a small plastic bag with water. Pectoral 
fin clips were taken and stored in 95% ethanol for DNA 
analyses. Once the fin clips were sampled, fish were allowed 
to recover in a bucket of well-oxygenated water from their 
habitat, and they were released at the location of capture. 
Genomic DNA was extracted using the DNeasy Tissue Kit 
(QIAGEN, Mississauga, ON, Canada), and DNA concen- 
tration was estimated using an ND-8000 NanoDrop spec- 
trophotometer (ThermoScientific, Ottawa, ON, Canada). 



DNA sequence data 

We used Polymerase chain reaction (PCR) to amplify two 
nuclear and one mitochondrial gene regions. The nuclear 
and mitochondrial loci have demonstrated utility in 
resolving question of population structure in other studies 
(e.g., Quattro and Jones 1999; Quattro et al. 2001). The 
nuclear loci are introns and although linkage to coding 
and/or regulatory variation is possible, they are likely to 
be neutral and any segregation of polymorphisms should 
reflect isolation of the populations not selection on these 
loci. Each 25 /.iL PCR reaction contained approximately 
0.06 ng of template DNA, 2.5 fiL of TopTaq Master Mix 
(QIAGEN), 2.5 /iL CoralLoad (QIAGEN), and 1.25 /iM, 
each, of specific forward and reverse primers. Reaction 
parameters were optimized for each primer set and are 
shown below. Reactions were checked for single amplifi- 
cation products (bands) using gel electrophoresis and 
ethidium bromide staining and sequenced at the Genome 
Quebec Innovation Center using an Applied Biosystems 
3730x1 sequencer (Burlington, ON, Canada). 

Creatine kinase 

The sixth intron in the Creatine Kinase- A (CK-A) nuclear 
locus was amplified using an oligonucleotide primer set 
previously optimized for various fish species, including 
cyprinids (Quattro and Jones 1999). The primers are as 
follows: CKA7F1: 5'-CCC AAG TTY GAG GAG ATC 
CTG AC-3' and CKA7R1: 5'-CCG TCG ACG ACC AGC 
TGC ACC TG-3'. These primers were designed to specifi- 
cally amplify the 6th intron from the CK-A locus, and 
not the CK-B locus, to eliminate the possibility of inad- 
vertently comparing paralogous sequences (Quattro and 
Jones 1999; Quattro et al. 2001). PCR conditions were: an 
initial denaturation of 2 min at 94°C followed by 35 
cycles of 1 min at 94°C, 40 sec at 48°C, 2 min at 72°C, 
and a final extension time of 7 min at 72°C. 

Triosephosphate isomerase 

The fourth intron of the Triosephosphate Isomeras A 
(TPI-A) nuclear locus was amplified using locus specific 
primers designed in a similar fashion to the CK primers (J. 
M. Quattro pers. comm.). Amplification of sufficient prod- 
uct to allow direct sequencing required semi-nested pri- 
mary and secondary PCR reactions. The primer sequences 
are as follows: TPIDF1: 5'-YTG ATY GGS CAG AAG GTG 
GC-3'; 302 - TPIB4 Rl: 5'-AGA ACC ACY TTR CTC CAG 
TC-3'; TPI4F2 - P301: 5'-GCA TYG GGG AGA AGC TRG 
AT-3'. The primary amplification was from TPIDF1 to 302 
- TPIB4 Rl. From this reaction, 2 /iL of product was 
diluted in 100 /.iL water and 2 fjL of this solution was used 
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as template in a second PCR reaction from TPI4F2 - P301 
to 302 - TPIB4 Rl. In both primary and secondary 
reactions, the PCR parameters were as follows: an initial 
denaturation of 2 min at 90°C followed by 44 cycles of 45 
sec at 90°C, 60 sec at 48°C, 60 sec at 72°C and a final 
extension time of 7 min at 72°C. 

16S mitochondrial locus 

A fragment of the 16S mitochondrial locus was amplified 
using "universal" primers: 16S AR: 5'-CGC CTG TTT 
AAC AAA AAC AT-3'; 16S BR: 5'-CCG GTC TGA ACT 
CAG ATC ACG T-3' (Simon et al. 1994). The PCR 
protocol was: an initial denaturation of 2 min at 90°C 
followed by 35 cycles of 1 min at 90°C, 1 min at 48°C, 
2 min at 72°C, and a final extension time of 7 min at 
72°C. 

Data analysis 

Raw sequence curve files were edited manually and 
aligned using CodonCode Aligner v 3.03 (CodonCode 
Corporation, Dedham, MA) with a final visual optimiza- 
tion of the alignment. The number of haplotypes (within 
and between location), haplotype diversity (Hd), nucleo- 
tide diversity (n), G ST and f ST _values were all calculated 
using DnaSP version 4.10 (Rozas 2009). Haplotype fre- 
quencies among all populations and distance matrices 
were calculated using ARLEQUIN version 3.11 (Excoffier 
et al. 2005). In order to determine allelic similarity and 
dissimilarity within and between populations, we used the 
Jost's D coefficient, as defined by Jost (2008) and calcu- 
lated in SPADE (Chao and Shen 2010). Significant dis- 
similarity was determined at a set value equal or above 
0.15 (Jost 2008). Maximum parsimony-based haplotype 
networks were estimated using TCS version 1.21 (Clement 
et al. 2000). Mantel Tests were conducted using XLSTAT 
(Addinsoft 2010) to correlate dissimilarity values with dif- 
ferences in dissolved oxygen values between locations. 

Results 

Creatine kinase 

We amplified and sequenced 198 bp of the sixth intron 
of the CK-A gene and identified 14 haplotypes across the 
seven collection locations (Fig. 2; Table 1) with an overall 
haplotype diversity of 0.791. (Table 1). The Njuguta River 
location had the highest diversity (Hd = 0.766), while the 
Dura River had the lowest diversity (Hd = 0.189). The 
maximum number of mutational steps between any of 
the haplotypes was four, although this maximum was 
only observed between a single individual from Mikana 
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Stream and the haplotype B. The most common haplo- 
type for this intron was haplotype A, found in 42 of 135 
individuals (31.11%, Fig. 2). This haplotype was observed 
in multiple individuals from all populations except the 
Dura River. In fact, all individuals from the Dura River 
site presented 1 of 2 unique haplotypes, C or D, not 
observed in any other population. 

Gst and Jost's D analyses suggest significant genetic dif- 
ferences between some of the populations at this locus 
(Table 2). F ST -values did not indicate significant diver- 
gence, but reflected and never contradicted G S T-values 
(data not shown). This pattern of P ST showing similar, but 
not significant, patterns of differentiation to those in the 
G S t analysis is consistent across all three loci and we have 
therefore presented only the G ST values. Values for Jost's D 
dissimilarity ranged from —0.112 to 1.000, indicating 
variability between the Dura River population and all other 
populations (Table 2). Furthermore, Jost's D values indi- 
cated that the Rwembaita Swamp population and the Inlet 
Stream West population were genetically different from 
each other with a dissimilarity value of 0.192. 

Overall, we found a significant positive relationship 
between geographical distance and both Gst and the 
Jost's D dissimilarity coefficient with populations that 
were farther apart being the most genetically differenti- 
ated (R 2 = 0.96, P < 0.0001 and R 2 = 0.97, P < 0.0001 
for Gst and Jost's D, respectively). Due to the large dis- 
tance between the Dura River population and all the 
other populations, additional analyses were conducted 
without the Dura River population. In the absence of the 
Dura River population, both G S t and Jost's D were asso- 
ciated with the geographical distance in pairwise compari- 
sons (R 2 = 0.33; P < 0.0001 and R 2 = 0.37; P < 0.0001 
for G S t and Jost's D, respectively). By contrast, we did 
not find a significant relationship between either genetic 
metric and dissolved oxygen at this locus, even when the 
Dura River was not included in the analysis. 

Triosephosphate isomerase 

We amplified and sequenced 150 bp from the fourth 
intron of the TPI-A gene and identified nine different 
haplotypes across the seven populations (Fig. 2; Table 1) 
with an overall haplotype diversity of 0.608 (Table 1). 
Similar to the results from the CK locus, the Rwembaita/ 
Njuguta Outflow Stream had the highest overall diversity 
(Hd = 0.763), while the Dura River was completely fixed 
at this locus with a haplotype diversity of 0.000. The 
maximum number of mutational steps between any of 
the nine haplotypes obtained was two. The most common 
haplotype for the TPI intron was haplotype B found in 
66 of 134 individuals (49%, Fig. 2) and in all locations, 
suggesting a single source population for this haplotype. 
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□ Dura River 

D Mikana Stream 

CI Rwembaita Swamp 

□ Inlet Stream West 
I Inlet Stream East 
H Njuguta River 

□ Rwembaita Outflow 
Stream 



Triosephophate Isomerase 



16S 





F# CD D 

Haplotype B was also the only haplotype observed in the 
Dura river population. The next most common haplotype 
was A, found in 55 of 134 individuals (41%) and in all 
populations except the Dura River population. 

Similar to CK, G S t and Jost's D analyses suggest signifi- 
cant genetic difference at this locus between some of the 
populations (Table 3). Values for Jost's D dissimilarity 
ranged from —0.065 to 0.794 and also showed significant 
differences between the Dura River population and all the 
other populations as well as between the Rwembaita 
Swamp population and all other populations except with 
the Rwembaita/Njuguta Outflow Stream population 
(Table 3). Additionally, the Njuguta River population and 
the Inlet Stream West population were marginally dissim- 
ilar from each other (Jost's D- value = 0.129). 

We observed a significant positive correlation between 
geographical distances and both genetic metrics 



Figure 2. Parsimony network of haplotypes 
for the Creatine Kinase intron (CK), the 
Triosephosphate isomerase Intron (TPI) and the 
16S gene. Coloured circles represent different 
haplotypes. Size of circles represent the 
number of individuals with a specific 
haplotype. Hollow nodes represent a single 
nucleotide polymorphism (SNP). Significance 
for connections is 0.95. 

(R 2 = 0.38, P < 0.0001 and R 2 = 0.46, P < 0.0001 for G ST 
and Jost's D, respectively). Nevertheless, when the Dura 
River was removed from the analysis, these relationships 
disappeared (P = 0.31 and P = 0.35 for G S t and Jost's D, 
respectively). By contrast, the relationships between the 
genetic metrics and dissolved oxygen between locations 
were significant when the Dura river was included 
(R 2 = 0.25, P = 0.04 and R 2 = 0.64, P < 0.0001 for G ST 
and Jost's D, respectively), but stronger when this popula- 
tion was excluded from the analysis (R 2 = 0.83, 
P < 0.0001 and R 2 = 0.83, P < 0.0001 for G ST and Jost's 
D, respectively; Fig. 3). 

16S mitochondrial locus 

We sequenced a 555 bp fragment of the 16S locus 
and identified seven haplotypes (Fig. 2) with an overall 
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Table 1. Sampling locations and genetic diversity indices for CK and TPI introns as well as 16S gene. 
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Env. Type, Type of environment; A/ Se q, Number of individuals sequenced; Poly, sites, Polymorphic sites; A/ Ha p, Number of haplotypes; Hd, Haplotype 
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fAII sites are normoxic except RS and Rwe, which are hypoxic. 



Table 2. G ST and Jost's D values for the CK intron. 
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0.354* 


0.376* 


Inlet Stream East 


1.000 




0.011 


-0.012 


0.003 


0.015 


0.026 


Inlet Stream West 


1.000 


0.073 




0.004 


0.008 


0.003 


0.041 


Mikana Stream 


1.000 


-0.061 


0.022 




-0.010 


-0.008 


-0.001 


Njuguta River 


1.000 


0.018 


0.050 


-0.0S4 




-0.021 


-0.003 


Rwembaita Outflow Stream 


1.000 


0.090 


0.020 


-0.038 


-0.122 




-0.004 


Rwembaita Swamp 


1.000 


0.092 


0.192 


-0.029 


-0.017 


-0.014 





G ST values are above the diagonal, the dissimilarity matrix is below the diagonal. Dura, Dura River; ISE, Inlet Stream East; ISW, Inlet Stream West; 
Mik, Mikana Stream; NR, Njuguta River; RS, Rwembaita Outflow Stream; Rwe, Rwembaita Swamp. 
*P< 0.05. 



haplotype diversity of 0.287 (Table 1). In contrast to the 
results from the nuclear loci, in this mitochondrial data 
set, the Dura River location had the highest haplotype 
diversity {Hd = 0.284). The Mikana Stream, Njuguta 
River, and the Rwembaita Swamp all had only a single 
16S haplotype and an Hd of 0.000. The most common 
16S haplotype was haplotype A found in 112 of 134 indi- 
viduals (83.6%). Haplotype B was observed in 17 of 134 
individuals (12.7%) and only in individuals from the 
Dura River. The remaining haplotypes were found only 
in single individuals. 

Similar to the two nuclear loci, G ST and Jost's D anal- 
yses suggest significant genetic differences at this locus 



between some of the populations (Table 4). In particu- 
lar, we observed a significant genetic isolation between 
individuals from the Dura River population and all 
other locations. Significant differences in G$r and Jost's 
D were also observed between the Mikana River and the 
Njuguta River, the Mikana River and the Rwembaita 
Swamp, and the Njuguta River and Rwembaita Swamp 
(Table 4). 

Similar to the results from the nuclear loci, there was a 
positive and strong correlation between geographical and 
genetic distance in the 16S data set (_R 2 = 0.36, P = 0.08 
and R 2 = 0.94, P < 0.0001 for both G ST and Jost's D, 
respectively). However, no correlation was observed 
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Table 3. G S j and Jost's D values for the TPI intron. 



Population 


Dura 


ISE 


ISW 


Mik 


NR 


RS 


Rwe 


Dura River 




0.228* 


0.184* 


0.188* 


0.283* 


0.323* 


0.624* 


Inlet Stream East 


0.278 




-0.017 


-0.019 


0.015 


-0.005 


0.094* 


Inlet Stream West 


0.177 


-0.045 




-0.025 


0.033 


0.017 


0.149* 


Mikana Stream 


0.202 


-0.049 


-0.065 




0.014 


0.005 


0.132* 


Njuguta River 


0.515 


-0.066 


0.129 


0.064 




-0.004 


0.081* 


Rwembaita Outflow Stream 


0.482 


-0.018 


0.050 


0.015 


-0.018 




0.033 


Rwembaita Swamp 


0.794 


0.210 


0.306 


0.301 


0.243 


0.075 





Gst values are above the diagonal, the dissimilarity matrix is below the diagonal. Dura, Dura River; ISE, Inlet Stream East; ISW, Inlet Stream West; 
Mik, Mikana Stream; NR, Njuguta River; RS, Rwembaita Outflow Stream; Rwe, Rwembaita Swamp. 
*P< 0.05. 
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Figure 3. Relationship between Jost TPI intron D values as function 
of the geographical distance along the waterway. Black points with 
circles represent the comparison between hypoxic and normoxic 
sampling sites; black squares represent comparisons between 
normoxic-normoxic sampling sites; black diamonds represent 
comparisons between hypoxic and hypoxic sampling sites. 

between geographical distance and either G ST (P = 0.363) 
or Jost's D (P = 0.834) when the Dura river was excluded 
from the analysis. Additionally, there was no correlation 
between the genetic metrics and dissolved oxygen even 
when the Dura River was removed from the analysis. 

Discussion 

Results from this study indicate that both geographical 
distance and dissolved oxygen regime contribute to 
genetic isolation and population genetic structure in 
B. neumayeri in the Kibale National Park, Uganda. Based 
on both nuclear and mitochondrial loci analyses, the 
Dura River population is isolated from the other 



sampling locations by geographical distance and landscape 
features, and shows the highest degree of genetic isolation 
in comparison with all the other locations. When the 
Dura River location was removed from the analyses, how- 
ever, the correlation between geographical distance and 
genetic differentiation disappears and a different pattern 
of differentiation is observed. 

Genetic structure of the B. neumayeri in the 
Kibale National Park 

Distance vs. oxygen regimes 

The relatively high level of genetic structure at the mito- 
chondrial locus between the Dura River and all other 
locations indicates little if any gene flow between the two 
drainages. Furthermore, the high dissimilarity between 
the Dura River and all other sample locations observed at 
both nuclear introns suggests that the isolation of the 
Dura River is not a relatively recent event. The clustering 
of the Dura River population in the haplotype analyses 
further supports the idea that there is little to no effective 
migration occurring between the Dura River and the 
other locations. Two isolated groups are readily apparent 
in the 16S haplotype network: the Dura River population, 
and a second group including all the other populations. 
The number of fixed mutations at the 16S locus is unli- 
kely to occur over a short period of isolation, suggesting 
that this isolation is not a seasonal event and one that has 
been established for many generations. 

Recent physiological studies (Martinez et al. 2004, 
2011) further support our hypothesis that genetic isola- 
tion of B. neumayeri is also associated with dissolved 
oxygen regime. Our results suggest reduced gene flow 
between populations from the two hypoxic Rwembaita 
locations and both the high oxygen locations and the 
normoxic Inlet Stream West populations, despite the fact 
that these two set of populations are only separated by 
328 m (one of the shortest distances between locations). 
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Table 4. G ST and Jost's D values for the 16S gene. 



Population 


Dura 


ISE 


ISW 


Mik 


NR 


RS 


Rwe 


Dura River 




0.658* 


0.664* 


0.740* 


0.729* 


0.664* 


0.723* 


Inlet Stream East 


0.941 




-0.013 


0.000 


0.000 


-0.013 


0.000 


Inlet Stream West 


0.941 


-0.003 




0.000 


0.000 


-0.013 


0.000 


Mikana River 


0.942 


0.003 


0.000 




1.000 


0.000 


1.000 


Njuguta River 


0.942 


0.003 


0.000 


0.000 




0.000 


1.000 


Rwembaita Outflow Stream 


0.941 


-0.003 


-0.006 


0.000 


0.000 




0.000 


Rwembaita Swamp 


0.942 


0.003 


0.000 


0.000 


0.000 


0.000 





G S t values are above the diagonal, the dissimilarity matrix is below the diagonal. Dura, Dura River; ISE, Inlet Stream East; ISW, Inlet Stream West; 
Mik, Mikana Stream; NR, Njuguta River; RS, Rwembaita Outflow Stream; Rwe, Rwembaita Swamp. 
*P< 0.05. 



This isolation suggests that habitat-specific selection 
pressures may contribute to reduced gene flow (or alter- 
natively, high gene flow may constrain local adaptation), 
as has been suggested in other studies (Hendry et al. 
2002; Crispo 2008). The TPI intron data further suggest 
that genetic differentiation is occurring between individ- 
uals from the hypoxic Rwembaita Swamp and individu- 
als from all the other sampling locations, except for 
individuals from the Rwembaita/Njuguta Outflow 
Stream, the only other hypoxic sampling location charac- 
terized by low- conductivity swamp water. Among all six 
sampling locations in the Rwembaita-Njuguta system, 
the Rwembaita/Njuguta Outflow Stream is the furthest 
from the Rwembaita Swamp (1997 m). In other words, 
the two hypoxic locations are genetically most similar 
while being geographically distant. The habitat in 
between these locations is a continuous papyrus swamp, 
possibly enabling swamp-dwelling fish to move more 
freely within the swamp system over broad distances. 
We observed strong correlations between DO and 
genetic dissimilarity at the TPI intron even when the 
Dura River is removed from the analysis, adding further 
support to our conclusion that oxygen regime is playing 
a significant role in determining genetic divergence and 
population structuring of B. neumayeri at micro- 
geographical scales. 

Interestingly, studies of the widespread cichlid P. multi- 
color in the Mpanga River system provided no evidence 
of population genetic structure with respect to oxygen 
regimes; however, there was evidence for isolation by dis- 
tance (Crispo and Chapman 2008). The different patterns 
of isolation seen in this study and the work from Crispo 
and Chapman (2008) may reflect differences in the spe- 
cies (i.e., cichlids vs. cyprinids), or differences in the level 
of resolution between the molecular loci (i.e., nuclear 
introns vs. microsatellites). Despite the fact that P. multi- 
color and B. neumayeri are from very phylogenetically 
distinct lineages, both occupy similar heterogeneous 
niches with respect to dissolved oxygen availability, and 



display similar patterns of divergence across DO gradients 
in morphological traits such as gill size. In P. multicolor, 
common-garden rearing studies have demonstrated a high 
degree of phenotypic plasticity in some traits related to 
oxygen uptake such as gill size, which may buffer the 
effects of local adaptation and allow for high rates of gene 
flow at the meta-population level (Chapman et al. 2008; 
Crispo 2008; Crispo and Chapman 2010). Further studies 
are needed to determine whether intron variation may 
exist across oxygen gradients in P. multicolor. While other 
environmental characters differ between low- and high- 
oxygen locations, several lines of evidence suggest that 
dissolved oxygen availability is driving divergence in phys- 
iological, behavioral, and morphological traits between 
B. neumayeri from low- and high- oxygen locations 
(Chapman and Liem 1995; Olowo and Chapman 1996; 
Chapman et al. 1999; Martinez et al. 2004, 2011). Our 
study therefore supports the hypothesis that non- 
geographical environmental factors can affect genetic vari- 
ability between subpopulations in a continuous system 
(Peichel et al. 2001; Hendry et al. 2002; Lecomte and 
Dodson 2004; Bekkevold et al. 2005; Funk et al. 2005; 
McCairns and Bernatchez 2008). In the case of B. neu- 
mayeri, it is likely that divergent selection has led to vari- 
ation in traits associated with oxygen uptake, in turn 
resulting in performance trade-offs for different respira- 
tory phenotypes that contributes to the maintenance of 
divergence in both phenotype and genotype (Schaack and 
Chapman 2003). The role of phenotypic plasticity in con- 
tributing to oxygen-related trait divergence in B. neumay- 
eri is an area for future study. 

G S t vs. Jost's D 

This study also allowed us to compare and contrast the 
results from distinct metrics of genetic structure. Some 
limitations that have been pointed out for Gst values 
(e.g., Charlesworth 1998; Nagylaki 1998; Jost 2008) are 
evident in our results. For example, in the analysis of the 
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16S data set, GsT-values revealed significant differences 
between B. neumayeri from the Rwembaita Swamp and 
both Mikana Stream and Njugata River; however, Jost's D 
showed no difference between these populations. It 
appeared that all but a single individual from the Dura 
River possessed a haplotype unique to this population, 
suggesting that G S i-values should approach 1 as opposed 
to 0.664-0.740 as indicated by our results. In comparison, 
Jost's D-values for the Dura River ranged from 0.941 to 
0.942, a better indication of strong genetic differentiation. 
Our results, then, appear to be consistent with concerns 
raised in earlier studies (Charlesworth 1998; Nagylaki 
1998; Jost 2008) that true genetic differentiation is not 
distinguishable when unique alleles exist in subpopula- 
tions, as well as when populations with different mean 
heterozygosities are compared. Similar differences were 
observed for the CK locus for which all the individuals of 
the Dura River presented one of four unique haplotypes 
indicating a very high differentiation (G ST -values = 0.339 
to 0.374 and Jost's D-values = 1.000). Such discrepancies 
may cause misleading conclusions with respect to popula- 
tion structure leading to inappropriate management 
implications and we therefore advocate for a better inves- 
tigation of the properties associated with different metrics 
used to estimate genetic differentiation. Furthermore, it is 
clear that G ST may yield less realistic measure of popula- 
tion differentiation than Jost's D and caution must be 
taken when interpreting and comparing the results from 
different studies. Recently, Raeymaekers et al. (2012) 
empirically and theoretically compared the performance 
of D and G ST in distinguishing population genetic structure 
and conclude that no single measure can quantify all 
aspects of genetic structure. The differences in resolution 
observed here between D and Gst are consistent with this 
conclusion. Interestingly, both their empirical and theoreti- 
cal work suggest that D is more useful in distinguishing 
longer term processes, and G ST in recent demographic 
events, while our empirical results appear to point in the 
opposite direction. This difference in conclusion suggests 
that either the separation of the B. neumayeri populations 
that we have studied is more complicated than we appreci- 
ate, or, that more work is necessary to completely under- 
stand the difference in performance of these two metrics. 

Conclusion 

Our study detected significant population differentiation 
in the African cyprinid B. neumayeri across dissolved oxy- 
gen regimes. We also observe strong genetic differentia- 
tion of the Dura River location from all other sites, likely 
a function of the large swimming distance and the pres- 
ence of geographical barriers between the Dura River 
location and all of the other samples sites. The remaining 
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six sampling locations, although geographically close to 
each other, present levels of genetic divergence associated 
with dissolved oxygen regime. Earlier studies on the bio- 
chemistry of this system found significant differences 
between different oxygen regimes (Martinez et al. 2004, 
2011) and demonstrated that these differences were not 
plastic, conclusions most easily explained by some degree 
of genetic isolation between populations. We find that 
within the Rwembaita Swamp system, the two hypoxic 
locations are genetically similar to each other, while the 
four normoxic locations are genetically similar to each 
other. Combined with the earlier biochemical studies, our 
results support the idea that aquatic dissolved oxygen can 
act as a selective force limiting gene flow among popula- 
tions of aquatic species and therefore should be consid- 
ered when implementing conservation plans and be of 
future concern when assessing environmental impact of 
human activities, especially on sensitive ecosystems such 
as tropical watersheds. 
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